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We ask whether the scaling exponents or the Kolmogorov constants depend on the 
anisotropy of the velocity fluctuations in a turbulent flow with no shear. According to 
our experiment, the answer is no for the Eulerian second-order transverse velocity struc- 
ture function. The experiment consisted of 32 loudspeaker- driven jets pointed toward the 
centre of a spherical chamber. We generated anisotropy by controlling the strengths of 
the jets. We found that the form of the anisotropy of the velocity fluctuations was the 
same as that in the strength of the jets. We then varied the anisotropy, as measured by 
the ratio of axial to radial root- mean- square (RMS) velocity fluctuations, between 0.6 
and 2.3. The Reynolds number was approximately constant at around i^A = 481. In a 
central volume with a radius of 50 mm, the turbulence was approximately homogeneous, 
axisymmetric, and had no shear and no mean flow. We observed that the scaling ex- 
ponent of the structure function was 0.70 =b 0.03, independent of the anisotropy and 
regardless of the direction in which we measured it. The Kolmogorov constant, (72, was 

also independent of direction and anisotropy to within the experimental error of 4%. 
f Email address for correspondence: gregory.bewley@ds.mpg.de 
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1. Introduction 

Theories of turbulence that give precise predictions apply by mathematical necessity 



to statistically isotropic flows. A conjecture, due to Kolmogorov (1941), is that flows 
characterised by large enough Reynolds numbers are locally isotropic, or isotropic at 
small scales, even in the presence of anisotropy at larger scales. Understanding the way 
turbulence tends to local isotropy from various states of anisotropy remains an important 
challenge in the study of turbulence. Its importance lies in the fact that almost all flows 
in the natural world and in technical applications are anisotropic, and not of sufficiently 
high Reynolds numbers for its influence to be neglected. Any useful theory or model of 
turbulence should therefore incorporate anisotropy, and the purpose of this investigation 
is to provide test cases for such theories. 



The conjecture of local isotropy has received much attention (e.g. Saddoughi & Veer- 



avalli|1994 Kurien fc Sreenivasan||2QQl Biferale fc Procaccia||2QQ5 ), and has been tested 



in flows with various kinds of anisotropy. Anisotropy can arise in different ways, for 



example through spatial variation of the mean flow (e.g. Tavoularis & Corrsin 1981), 
or through anisotropy in the velocity fluctuations. These different forms of anisotropy 
might have different effects on the small scales of turbulence. No test has been performed 
previously to isolate the influence of anisotropy in the fluctuations only, without also 
introducing shear. Our task was to create a device with which we could do this. We show 
that it is possible in a single apparatus to produce turbulence with negligible shear, and 
with a chosen level of anisotropy on the large scale. With the apparatus, we tested the 
validity of the hypothesis of local isotropy under increasingly anisotropic conditions. 

We produced anisotropy in the velocity fluctuations by introducing an asymmetry in 
the agitation of the turbulence. Evidence suggests that there is a causal relationship be- 
tween asymmetry in the agitation of turbulence, and anisotropy in the energy-containing 
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scales of turbulence. For example, asymmetry and anisotropy are present at the same 



time in turbulent jets (e.g. Hussein, Capp fc George|1994 ), Taylor- Couette flows (e.g. An- 



[dereck, Liu fc Swinney||1986 ), and wind tunnels with specially designed shear generators 
(e.g. Tavoularis fc Corrsin||198T ). Turbulence produced by computer simulation can also 



be forced in an asymmetric way; Yeung fc Brasseur (1991) studied the influence of this 
asymmetry. One exception is the turbulence produced by a grid in a wind tunnel with 
a contraction, where the axis of the tunnel introduces an asymmetry, yet the turbulence 
produced can be nearly isotropic, but decaying ( Comte-Bellot fc Corrsin||1966 ). There is 
also a pattern in the relationship between asymmetry of the forcing and anisotropy of 
the turbulence. Machines with a single axis of symmetry, such as the von Karman flow 



(e.g. Ouellette et al. 2006) and wind tunnels (e.g. Staicu, Vorselaars fc van de Water 



2003), produced axisymmetric turbulence. Machines with more axes, such as the one 



developed by Hwang fc Eaton (2004), produced isotropic turbulence. A careful study by 



Zimmermann et al. (2010) showed that six axes were sufficient to produce turbulence 



without a preferred direction. Our machine had 16 axes. 

We generated both isotropic and axisymmetric flows. The apparatus worked by mod- 
ulating the relative intensity of 32 mixers distributed over a sphere. When all the mixers 
were driven with the same intensity, isotropic turbulence resulted. We introduced a pre- 
ferred axis by driving mixers near one of the 16 axes of the machine either more strongly 
or more weakly than the other mixers. This produced anisotropy in the turbulence. The 



geometric form of the apparatus was inspired by the one of Zimmermann et al. (2010), 
except that ours used 32 mixers instead of their 12, and air as the working fluid instead 
of water. Our mixers were loudspeakers coupled to nozzles; the diaphragm of the loud- 
speaker pushed and pulled air through an orifice, which formed a turbulent jet. This 



technique for mixing air was invented by Hwang fc Eaton (2004), though our imple- 
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mentation differs in several details, described below. We used laser Doppler velocimetry 
(LDV) to measure the fluid velocity near the centre of the apparatus. 

We measured the influence of anisotropy on velocity structure functions, which are the 
moments of velocity differences across two points separated in space. The second order 
structure function is 

Ai(r-) = {{Ui{r) - Wi(0))(w,(r) - u,(0))) , (1.1) 

where Ui{r) is the fluctuation of the z-th component of the velocity at location r and 
(•) denotes temporal averaging. Local isotropy demands that for small enough r and 
large enough Reynolds number, the structure function is isotropic. In this case it can be 
written in terms of a single independent scalar function. 



Kolmogorov (1941) surmised that in the inertial subrange, where he thought that 
neither the viscosity nor the large-scales would influence the dynamics, the structure 
functions follow a power law as a function of spatial separation. For representative parts 
of the structure function tensor given in equation |1.1[ this can be written 

DUri) ^ {{uM) - uM?) = \ C^2^ i^r^f'^' , (1-2) 

Dr,,,{z) ^ {{UrM-UrM?) = ^0^'^ z)^'-'' . (1.3) 

These functions are called transverse because the velocities, and , are perpendicular 
to the vectors pointing from to ri and from to respectively. We have introduced a 
certain coordinate system that we describe later. Here, e is the energy dissipation rate per 

(r) (z) 

unit mass, the dimensionless constants, C2 and C2 , are called Kolmogorov constants. 



and the scaling exponents, ^2^^ and ("2^^ both equal 2/3^^^ according to Kolmogorov 



( 1941 ) for locally isotropic turbulence. The Kolmogorov constants are also equal to each 
other if the turbulence is locally isotropic, but their value must be determined by ex- 
periment. It is interesting to note that the corresponding Kolmogorov constants for the 
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Lagrangian structure function vary with angle by more than 20% at a Reynolds numbers 



of 815 in an anisotropic von Karman swirling flow (Ouellette et al. 2006). Our work 
addresses the angular dependence of the Eulerian C2, and is in part an extension of the 
above work. 

One of our measures of anisotropy is the ratio of Kolmogorov constants measured 
in the two different directions, 02^"^ /C2^\ Variation of this quantity could account for 
differences between published values of the Kolmogorov constant. Such differences are 



evident in Sreenivasan (1995), for example, where values for the constant were collected 
from various experiments and found differ from each other with a standard deviation 
of about 10%. Could it be that different experiments measured different values for the 
Kolmogorov constant because they had different levels of anisotropy, or else because they 
measured the constant in different directions? In our study, we found no dependence of 
the Kolmogorov constants on the anisotropy, as shown in section |4] 

We also report on the difference between the scaling exponents measured in the two 
directions, A(2 = (2^^ ~ Co!^^ - Here, we sought to interpret the observation made by 



Shen & Warhaft ( 2002 ) that in anisotropic wind tunnel flows the scaling exponent in the 



cross-tunnel direction was about 0.1 smaller than the exponent in the direction of the 
mean flow. In our shear less, axisymmetric flow we found, in contrast, that the difference 
between the scaling exponents, AC2, was zero to within the measurement accuracy. As 
is explained in section [4j we measured the exponents using the Extended Self-Similarity 



(ESS) technique introduced by Benzi et al. (1993). 

In the following section [2j we describe the apparatus and techniques. In section [3) we 
show that we could control systematically the anisotropy of the turbulence, and report 
on the extent to which the turbulence had no shear and was axisymmetric. A reader 
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(a) (b) 

Figure 1. (a) The jet generator, (b) The layout of the turbulence chamber and the LDV probes 
(black cylinders), one near the bottom of the image, and one to the right. The one aligned with 
the xi-axis measures the components of particle velocities in both the X2 and xs directions. The 
probe aligned with the X2 axis measures velocities in either the xi or the X3 direction. 

with no interest in experimental detail can proceed to section [4j where we report on the 
measurements of structure functions. We conclude in section HI 

2. The Experiment 

As shown in figure [l] (b), the turbulence chamber had the shape of a truncated icosa- 
hedron and a diameter of 99 cm. It was made of wood 

joined together with nylon straps and glue. In the centre of each face a circular hole 
was cut for the jet generator. In addition, further holes adjacent to the jet generators 
were cut for optical access. As the shape is that of a soccer ball we will use this term 
from now on to describe the turbulence chamber. 

As depicted in figure Ilia), the mixers were jets, similar in principle to those developed 
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for flow control applications (e.g. Glezer & Amitay 2002). A 150 W loudspeaker with 
a diameter of 16.5 cm and a uniform frequency response between 50 and 3000 Hz was 
mounted on each face of the soccer ball and pointed towards its centre. The loudspeakers 
were coupled to conical nozzles with opening angle 30°, length 4.3 cm, and orifice diameter 
5 cm. When backing away from the orifice the diaphragm drew in air from all directions, 
while it blew it out in a single direction when moving toward the orifice. By driving the 
speaker sinusoidally, we generated a pulsating turbulent jet. The sound level inside the 
soccer ball was typically 135 dB, which corresponded to 0.2% of the measured turbulent 
kinetic energy. Thus we expect the sound to have an negligible effect on the turbulence 
and the measurements of velocity. 

The fiow chamber was similar to the ones described by [Hwang fc Eato"n (2004), Web- 



ster, Brathwaite & Yen (2004) and Lu et al. (2008), who employed cubic Plexiglas boxes 



from which the corners were cut off. A speaker-driven jet was mounted at each of the 
eight cut-off corners, and pointed toward the centre of the chamber. The fiow chamber of 



Warnaars, Hondzo & Carper (2006) was a rectangular Plexiglas box with two speakers 



placed at each square end of the box. A grid was placed in front of each speaker to excite 
small-scale turbulence. We compare the fiow parameters in these apparatuses, as well as 
the one of |Goepfert et a/.| ( |2010D , in table [l] 

We chose the signals that drove the loudspeakers to generate the desired anisotropy 
and low mean velocity in the fiow. Figure [2] shows the power spectrum of the signal. The 
50 Hz amplitude-modulated sine signals were digitally generated at a sampling frequency 
of 3 kHz. All speakers were driven in phase, and 50 Hz produced the strongest jets. The 
50 Hz oscillation was not, however, detectable in the spectrum of the fiuid velocities 
measured near the middle of the soccer ball. We modulated the amplitude of each driving 
signal with noise ( Fox et a/. 1 1988 ), though other forms of noise would probably also work. 
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Hwang 

& Eaton (2004) 


Webster 
et al (2004) 


Warnaars 
et al. (2006) 


Lu et al. 
(2008) 


Goepfert 
et al. (2010) 


Present 
work 


Geometry 


Truncated 
cube 


Truncated 
cube 


Rectangular 


Truncated 
cube 


Octahedron 


Truncated 
icosahedron 


Flow medium 


Air 


Saltwater 


Water 


Air 


Air 


Air 


a (ms~^) 


0.87 


0.0097 


0.00069 


0.60 


0.88 


1.1 


O'u/O'v 


1.03 


1.00 


0.98 




0.95 


0.94 


U/(7 


0.022 


0.07 






0.03 


0.04 


e (m^s-^) 


11 


2.5 X 10-^ 


1.25 X 10-^ 


1.4 


5.8 


6.7 


77 (^im) 


130 


450 


950 


220 


155 


155 


Trj (ms) 


1.2 


200 


1000 


3.3 


1.6 


1.5 


Rx 


220 


68 


5.9 


260 


250 


481 



Table 1. Chamber geometry and flow statistics from various speaker-driven flows. Symbols 
from top to bottom: RMS velocity fluctuations (a), ratio of RMS velocity fluctuations {(Tu/cfv)^ 
ratio of mean flow to RMS fluctuations (U/a), energy dissipation rate (e), Kolmogorov length 
scale (77), Kolmogorov time scale (rry), Taylor Reynolds number (Rx). For each study, we give 
the highest Reynolds number, lowest mean flow and best isotropy reported by the authors. For 
our apparatus, we report only on data acquired under isotropic forcing. 

The correlation time was 0.1 seconds, which was approximately equal to the large-scale 
eddy turn-over time, L/cr, L being a characteristic length scale describing the large-scale 
motions of the flow, and a being the root mean squared (RMS) velocity fluctuations. This 
condition ensured that fluctuations in the energy input rate to the turbulence occurred on 
time scales equal to or faster than the turbulence decay time, so that the turbulence was 
in a steady state. In any case, we found that the statistical properties of the turbulence 
were insensitive to the correlation time. 

To avoid mean flows, the amplitudes of each driving signal were adjusted so that the 
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>-4 




Figure 2. The power spectrum of the voltage apphed to the loudspeakers. 



sum of the amplitudes of all signals was zero. This also reduced the amplitude of the 
sound generated and the air exchange between the inside of the ball and the room. 

We selected the desired anisotropy by separately adjusting the RMS amplitude of each 
speaker As we restricted ourselves to cylindrically symmetric driving, the asymmetry 
could be characterised by the ratio of the axial amplitude, aaxiai, to the radial amplitude. 



A forcing is then described as exhibiting oblate spheroidal asymmetry when < A < 1, 
spherical symmetry when A = 1 and prolate spheroidal asymmetry when A > 1. 

As we have 32 speakers distributed with icosahedral symmetry, the RMS amplitude for 
a given loudspeaker was set according to the following scheme. We defined an ellipsoid 
that had the same centre as the soccer ball, and whose major axis to minor axis ratio was 
A. We then calculated the intersection between the surface of the ellipsoid and the vector 
that pointed from the centre of the soccer ball to the centre of the given loudspeaker. 
The distance between this intersection and the centre of the soccer ball set the relative 
RMS amplitude for the given loudspeaker. 



^radial* 



A = 



^^-axial 



(2.1) 



^radial 
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We characterised the flow using a TSI laser Doppler velocimetry (LDV) system. We 

used oil droplets as tracer particles, generated by a Palas AGF 10.0 aerosol generator. 

The droplets had a most probable diameter of about 3 |Ltm, and a volume fraction of less 

than 10~^. Droplets of this diameter settle in still air at 200 iitm s~^, which was much 

smaller than their RMS velocities in the turbulent flow. These oil droplets were passive 

tracers, since the Stokes number in our flow was very small (approximately 0.02) (e.g. 



Bewley, Sreenivasan & Lathrop 2008) and the accelerations of the turbulent flow were 



much larger than the acceleration of gravity (e.g. Voth et ar||2002 ). 

As shown in figure [TJb), we used two LDV probes. One probe measured two compo- 
nents, and the other only one component of the velocity of individual tracer particles. 
The measurement volume was approximately ellipsoidal in shape, and was approximately 
100 |LLm in diameter and 2 mm in length, which allowed us to resolve scales larger than 
the dissipation scale of the turbulence. The mean data sampling rate was between 300 
and 3000 samples per second, depending on the probe. We typically collected 2 x 10^ 
data points per aspect ratio of the forcing and per spatial position of the probes. 

LDV data are known to suffer from biases, since the system observes particles with 



larger velocities more often than those with smaller velocities (see e.g. Albrecht et al. 



2003). We used the residence time weighting suggested by Buchhave (1975) and Buch- 



have, George Jr. & Lumley ( 1979) to correct our single point statistics. We did not correct 



for possible bias in our two-point statistics. To calculate the two-point statistics the ve- 
locity signals were first processed using a slotting technique, i.e., we searched for samples 
coming from each of the two probes with time separations falling within a temporal bin 
of 1 ms duration, and within each temporal bin, the data from each probe was averaged. 
Approximately 10^ — 10^ points contributed to the statistics for each spatial separation. 
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We found that as long as the temporal bin was smaller than 15 ms the difference in the 

results was less than 1%. 

Before we proceed to the measurements, we describe the coordinate system, which was 
fixed to the axis of symmetry of the turbulence. Please note that the coordinate system 
was not fixed in the laboratory frame, as is shown in figures [l] and [4j Axisymmetric 
turbulence has two principal axes. However, our measurement apparatus measured ve- 
locities only at points that lay along a single line in the laboratory frame. In order to 
sample the statistics of turbulence along both axes, we rotated the axis of symmetry of 
the turbulence by taking advantage of the symmetries of the soccer ball. 

The soccer ball was oriented such that one of the symmetry axes of the forcing lay 
along = 107° and (j) = —4°, which is labelled as x[ in figure [4j and is close to the xi- 
axis. The other symmetry axis lay along ^ = 4° and (j) = 180°, which is labelled as x'^ in 
figure[4j and is close to the xs-axis. Because the primed and unprimed coordinate systems 
are close to each other, we do not distinguish between them in the rest of this paper. 
This simplification does not impact our conclusions, since we find no difference in the 
small-scale statistics between the two directions, despite the fact that they approximately 
orthogonal. 

Figure [3] shows the coordinate system, (ri, r2, 2:), which was aligned with the symmetry 
axis of the forcing. The coordinate system has two orientations with respect to the 
laboratory coordinate system, corresponding to the two cases described above. Hereafter, 
'axial' refers to both the direction of the velocity component measured along the axis of 
symmetry, and, in discussing two-point statistics, separations along the axis of symmetry. 
Similarly, 'radial' refers to both the direction of the velocity components perpendicular 
to the axis of symmetry, and to separations perpendicular to the axis of symmetry. In 
addition, we refer to two-point quantities whose separation vector lies along the axis of 
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z xs 



uzin) 




Ur2{z) 



Wr2(0) 



r2 ^3 




Ur2{ri 





Uri{z) 




r2 X2 



n X2 



(a) 



(b) 



Figure 3. Orientations of the body coordinate system of the forcing with respect to the 
laboratory frame, and the conventions we use (a) when the axis of symmetry of the forcing lay 
close to the xs axis of the laboratory frame, and (b) when the axis of symmetry lay close to xi 
axis of the laboratory frame. The coordinate system (xi, X2, xs) is fixed in the laboratory frame, 
and the coordinate system (ri,r2, z) is fixed with respect to the symmetry of the forcing. 

symmetry as axial quantities, and those with radial separations as radial quantities. For 
example, fixi) is denoted as f{z) when the axis of symmetry of the forcing lay along 
xi, and is called an axial quantity. In keeping with the geometry described above, our 
'axial' measurements were in reality about 15° away from the axis of symmetry, and the 
'radial' measurements were about 4° from its normal. 

In the first part of this paper, we discuss measurements made at a single point of all 
three orthogonal components of velocity. For these measurements, both probes observed 
a fixed volume (aligned to within 10 iitm) close to the centre of the soccer ball, which is 
the point (0, 0, 0) in the coordinate system shown in figure[4j One probe measured a single 
component of the particle velocities, namely (0, 0, 0). The second probe measured two 
orthogonal components of the particle velocities, namely Ux2{^-, 0, 0) and Uxr^{^^ 0, 0). For 
measurements described in the second part of this paper, we aligned the probes in a 
similar way, except that we positioned the second probe at different stations along the 
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Figure 4. The schematic shows the coordinate system in the laboratory frame, xs was aUgned 
with the vertical. The axis of symmetry of the forcing lay along either x[ or X3, depending on 
the experiment. Our measurements were made at points along xi. 

Xl axis, using a programmable linear traverse. This probe now measured Ux2{xi,0,0) 
and Uxq{xi^O^O). In addition, the single- velocity-component probe was rotated 90° to 
measure iia^3(0, 0, 0), which was coincident with one of the components measured by the 
two component probe when Xi equaled zero. 



For each of a series of forcing anisotropics. A, given by equation 2.1, we collected two 
data sets. For one data set, we aligned the axis of symmetry of the forcing close to X3, 
and for the other, we aligned the symmetry axis close to xi. For each value of the large- 
scale anisotropy, and each orientation of the symmetry axis, we collected data with the 
two-component LDV probe stationed at various positions along the xi-axis. We varied 
the anisotropy while fixing the quantity 
1 



K = 



{ui^ (0, 0,0)) + « (0, 0, 0)) + « (0, 0, 0)) 



(2.2) 



where the fluctuating parts of the velocities, are denoted by and (•) denotes 
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Axis Orientation 

Velocity Statistics (m s~"^) 

z X 



t/xi (0,0,0) 


0.099 


0.14 


[7.2(0,0,0) 


0.0055 


-0.0044 


[7.3(0,0,0) 


-0.038 


0.027 


f^xi (0,0,0) 


0.99 


0.99 


(0,0,0) 


0.98 


0.98 


(0,0,0) 


0.96 


0.98 



Table 2. Three- component velocity statistics for symmetric forcing, measured at the centre 
of the soccer ball. The symbols U indicate the mean flow and the symbols a indicate the RMS 
velocity fluctuations. 



temporal averaging. In other words, we fixed the turbulence kinetic energy in the centre 
of the ball. 



3. Anisotropy, axisymmetry and homogeneity 

Three-dimensional velocity measurements made at single points in the centre region 
of the soccer ball show that the flow was approximately homogeneous and axisymmetric. 
Furthermore, the anisotropy of the turbulent velocity RMS fluctuations followed the 
anisotropy of the forcing signal. For the case of spherically symmetric forcing, with A = 1 



(see equation 2.1), we expected and indeed found that the turbulence was isotropic. In 
this case, the components of the fluctuating velocity were each approximately 1 m s~^, 
and the mean velocity was less than 0.15 m s~^, or less than 15% of the fluctuations, as 
shown in Table [2j 

Figure |5] shows the ratio of axial to radial velocity fluctuations as a function of the 
forcing anisotropy, A. Hereafter, we refer to ('uf (^, ri, r2))"^/^ as cr^(z, ri, r2), and cr^ as 
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Figure 5. Velocity fluctuation ratios as a function of the forcing anisotropy, A. O is az/^n 
and A is Gzl^r-i- Data were collected with the axis of symmetry lying along X3. We obtained 
similar results (not shown) when the axis of symmetry lay along x\ . 

(7^(0,0,0). Clearly the anisotropy of the turbulence followed the variation of A. In addi- 
tion, the anisotropy was nearly the same whether measured in the r\ or r2 directions. 
Thus the turbulence was close to cylindrically symmetric, though this was less so at 
extreme values of A. The degradation in the cylindrical symmetry might be explained 
by the fact that the number of loudspeakers doing most of the work to drive the tur- 
bulence decreased as the value of A moved away from one. This is the nature of the 
forcing algorithm described before. As the number of loudspeakers effectively decreased, 
the turbulence probably became more sensitive to mechanical differences between the 
speakers, and to misalignments of the nozzles. 

Next, we consider the degree to which the axisymmetry was spatially uniform. To test 
this, we evaluate how the velocity fluctuation ratios vary with distance from the centre 
of the soccer ball. Figure |6|a) shows the ratio of the two radial fluctuating velocities, 
(7^2(0, 0, zj/dri (0, 0, 2:), at various distances from the centre of the ball along the axial 
direction. Within 50 mm, the values of this ratio deviated by less than 10% from 1, 
which indicates that the turbulence was close to cylindrically symmetric. Figure IgF 
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A 


0.33 


0.4 


0.5 


0.8 


1.0 


1.25 


2.0 


2.5 


(Jz (ni s""*") 


0.74 


0.75 


0.79 


0.91 


1.08 


1.21 


1.42 


1.47 




1.26 


1.24 


1.27 


1.19 


1.15 


1.03 


0.83 


0.75 


G z / CTro 
^ / '2 


0.59 


0.60 


0.63 


0.77 


0.94 


1.17 


1.71 


1.97 


A (m s ) 


1.87 


1 o o 


1 no 




1.91 


1.79 


1.d9 


1 (id 


\Uz/Ur,\ 


0.31 


0.16 


0.13 


0.15 


4.04 


2.57 


7.58 


23.9 


\U\ (m s-^) 


0.17 


0.18 


0.13 


0.10 


0.05 


0.07 


0.14 


0.24 


{Uri Ur2)/ (C^ri CT^s ) 


-0.07 


-0.06 


-0.06 


-0.07 


-0.07 


-0.05 


0.01 


0.04 


Ur2)/{<yz CTrs) 


-0.06 


-0.07 


-0.06 


-0.07 


-0.07 


-0.06 


-0.01 


0.02 




5.4 


4.9 


5.4 


6.0 


6.7 


7.0 


6.5 


6.4 


ry (^im) 


163 


167 


163 


159 


155 


153 


156 


156 


Tr^ (ms) 


1.7 


1.8 


1.7 


1.6 


1.5 


1.5 


1.6 


1.6 



Table 3. The table shows the turbulence parameters, for various loudspeaker RMS ampli- 
tude ratios (A). The quantity K is the turbulent kinetic energy, ^(cr^ + 2(7^2)5 1^1 is 
{U^ + 2U^^y^^. Data were collected at the centre of the soccer ball when the axis of symmetry 
was horizontal. We obtained similar results (not shown) when the axis of symmetry was vertical. 

shows the ratio of the axial fluctuating velocity to one of the radial fluctuating velocities, 
^z{^i^ 0)/c'"r2(^i7 0)- The values at ri = correspond to the values shown in figure [sj 
Again, the velocity fluctuation ratio is approximately constant within 50 mm of the 
centre of the ball. We infer that the anisotropy was approximately uniform within a 
central region with radius 50 mm. 

To evaluate the homogeneity of the turbulent fluctuations, we show measurements of 
(7^2(0, 0, z)/crr2(0, 0, 0) and az{ri,0,0)/az{0,0,0) in figures ^a) and (b). These quan- 
tities compare the amplitudes of the velocity fluctuations at locations away from the 
centre to those in the centre. Within a radius of 50 mm, the difference between 
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Figure 6. To demonstrate homogeneity of the turbulence axisymmetry, we show (a) the ratio of 
the two radial components of the velocity fluctuations, (0, 0, z)/ari (0, 0, z), and (b) the ratio 
of the axial fluctuations to the radial fluctuations, cr;z(ri, 0, 0)/crr2 (^i, 0, 0), at various stations 
moving away from the centre of the ball. The values in the legends, and in all other legends in 
this paper, are those of the data in (b) for ri = 0. 
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Figure 7. To demonstrate homogeneity of the turbulent fluctuations, we show (a) the radial 
fluctuations at positions along the axis, (7^2 (O5 O5 ^)/^r2 (O5 O5 0)? ^-nd (b) the axial fluctuations at 
positions normal to the axis, az{ri,0,0)/az{0,0,0). As in figure [g] the values in the legend are 
the anisotropics measured at the centre of the soccer ball. 

(7^2(0, 0, z)/crr2(0, 0, 0) and its value at the origin was within 5% of the value at the 
origin. The inequality also held for ^^(ri, 0, 0)/<j2(0, 0, 0). 
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Figure 8. The mean flow in two directions as a fraction of tfie fluctuating velocities in those 
directions. The values in the legend are the anisotropies measured at the centre of the soccer 
bah. 

Anisotropy can manifest itself not only in the fluctuations, but also through spatial 
variation of the mean velocity. In figures|8|a) and (b), we show the ratios of the mean flow 
velocities to the velocity fluctuations, (0, 0, z)/ (Jr^ (0, 0, z) and Uz (ri , 0, 0)/ (ri ,0,0). 
The shape of the curves are difficult to interpret, but may indicate that the natural centre 
of the turbulence was offset from the centre of our coordinate system. We concluded that 
variation in the mean velocity was negligible because it was typically more than ten times 
smaller in magnitude than the fluctuations, so that the energy in the mean flow was less 
than 1% of that in the fluctuations. Only at the extreme values of the anisotropy was 
the mean velocity as much as 15% of the fluctuations, which may have been due to the 
sensitivity of the turbulence to small differences between the speakers in this range of 
anisotropies, as discussed above. 

Turbulence production by a mean shear is gauged by the Reynolds stresses, or the cross 
correlations between orthogonal velocity components. We measured two such stresses, 
{ur^ (0, 0, z) Ur^{{)^ 0, z)) and (^/^(^i, 0, 0) Ur^{ri^{)^ 0)), at various locations, z and ri, along 
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the two axes. The stresses were less than 7% of the corresponding kinetic energies, (7^ cr^s 

and az cf^-^ , respectively, and so had a only a small influence on the flow. 



4. Universality in the second order velocity structure functions 



We first use the transverse structure functions, defined by equations L2 and |1.3[ to 
estimate the Reynolds number of the turbulence. In figure |9) the structure functions 
Dzzi^i) and Dr2r2{^) normalized by the Kolmogorov scaling. That is, we solve the 
equations for e, the dissipation rate. Here, we assume that the Kolmogorov constant 



has the value given by Sreenivasan (1995), C2 = 2.1, which is the mean of values taken 
from a collection of experimental studies. We assess the validity of this assumption later. 
We take the peak value of the function as our definition of the dissipation rate of the 
flow. The values estimated from 1)^2^2 (^) given in table [s] The corresponding values 
estimated from Dzzi^i) differ from those estimated from Dr2r2{^) by no more than 13%. 
The dissipation rates depended slightly on the driving and reached a value of 6.7 m^s~^ 
for the case of isotropic turbulence. Under isotropic forcing, the Taylor scale was then 



mm, where u is the kinematic viscosity of air at room temperature 
(1.57 X 10~^ m^ s~^), and the spherically averaged RMS velocity fluctuations was then 



a = y(^^+^a^)73, which gave a Taylor-microscale Reynolds number of Rx = Xa/u = 
481. The Reynolds numbers for the anisotropic cases were similar. 

We defined a scale-dependent measure of anisotropy, Dzz{ri)/Dr2r2{^) where ri = z. 



The data are shown in figure 10 Local isotropy requires that Dzz/Dr2r2 approaches 
one at small scales. In the limit of large separations, the ratio should approach cr^/cr^, 
since the velocities Ui{r) and Ui{0) are uncorrelated when r is large enough. Although 
we did not observe this regime, the values of the ratio in anisotropic cases did separate 
from the isotropic value as if to approach the values cr^/a^ at larger, unobserved, scales. 
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Figure 9. In (a) and (b), the axial structure functions, Dr2r2 {z)^ cind radial structure functions, 



Dzz{ri), are normalized by Kolmogorov's inertial range predictions given by equations 1.3 and 
1.2 In order from bottom to top, the anisotropy, az/cTr, was 0.59, 0.77, 0.94, 1.16 and 1.98. 
Each curve was shifted by 1.5 m^s~^ with respect to the one below it (except the bottom one). 

The figure shows that we observed local isotropy, within the experimental error, for 
separations smaller than about 50 mm. This approach to isotropy, then, is a consequence 
of the velocity correlations inherent in the dynamics of turbulence. We extract further 
information from the structure function ratio through two interpretations, described 
below. 

Because the structure functions do not exhibit clear scaling, as is evident in figure |9| 
we employ ESS to extract scaling exponents. We employed it in a new way, to compare 
the structure functions measured in different directions, rather than to compare structure 
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ri , z (mm) 

Figure 10. The ratios between structure functions in different directions but in the same flow, 
Dzz(ri) / Dr2r2{z)i approach the isotropic value of one as the separation n — z decreases. Note 
that noise became important below separations of about 4 mm, as can also be seen in figure [9] 
Results for data taken at other values of the large-scale anisotropy are consistent with those 
shown here. The values in the legends are the anisotropy, GzjcFr^ measured at the centre of the 
soccer ball. 



functions of different orders. According to Benzi et al. (1995), ESS implies that 



(j) 



(4.1) 



where L is a large scale of the flow and A'^^^ are unknown constants. For standard ESS, 
one finds that the function f{xi/r]) is independent of the order of the structure function. 
We assumed that it is also independent of the direction in which the structure function is 
measured. In this way, we assumed isotropy of the function /, but not of the underlying 
scaling. It follows that 



(z) — ^ ^1 — ^ ^1 



(4.2) 



whenever ri equals z. The ratio is a form of ESS because it eliminates the influence of 
the function / in order to reveal the underlying scaling set by the exponents ("2"^^ To 



determine we fit power laws to the data in figure 10 in the region between 4 and 
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Figure 11. In (a), the O symbols are the scahng exponent differences measured through 
the use of equation |4.2| For the □ symbols, the differences were computed through the use of 
equation |4.3| Two data points above the top of the plot and are out of view in order to make 
the deviations of the other points more visible. The values of these two points are 0.5±0.4 and 
0.3±0.3 in order of increasing anisotropy. The exponents measured using equation |4.3| are shown 
in (b). Here, the O and □ symbols mark radial and axial scaling exponents, respectively. 

40 mm, the lower bound being set by the appearance of noise in the data, and the upper 
bound by the emergence of the large-scale cutoff. It can be seen in figure pTj^a) that the 
dependence of AC2 on the anisotropy is non-monotonic. 



Figure 11 ^a) also shows exponent differences derived from standard ESS. In the original 
expression of ESS, 



(4.3) 



where Dm is the third moment of the absolute value of the velocity difference, {\ui{xj) — 
Ui{0)\^). For ESS to be present. Dm = Bi (r/L) f{r/r]). In figure 12, we plot D^z against 



(r) (z) 

Dzzz^ and Dr2r2 against Dr2r2r2^ order to uncover the scaling exponents Q and Q , 



respectively. We fit power laws to the data in figure 12 for r and z between 4 and 40 mm, 
as before. The exponents of these power laws are shown in figure [njb). Note that for 
two exponents measured in the radial direction of prolate spheroidal turbulence, the 
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Figure 12. In (a) and (b), the axial structure functions, Dr2r2{z), and radial structure func- 
tions, Dzz{ri), are plotted against their third-order counterparts. For each curve in order from 
bottom to top, the anisotropy, az/cTr, was 0.59, 0.77, 0.94, 1.16 and 1.98. Each curve was shifted 
by a factor of 2 with respect to the one below it (except the bottom one). The solid lines are 
power laws with an exponent of 0.70. The uppermost curve in (b) is marked with stars in order 
to make the data distinguishable from the other data. For these data, the third order structure 
function did not converge. 

data were insufficient to bring about convergence of the third order structure functions. 



This is evident in the scatter present in the corresponding data in figure 12 We note 
speculatively that for these two conditions, the Reynolds stresses were about 0.08 higher 
than for all of the other cases, as seen in table [3| and that this may point to a role for 
shear in altering the scaling. 



As seen in figure ll[b), the dependence of the scaling exponents, ^2^^ and ("2^^ on 



anisotropy was non-monotonic. Furthermore, deviations of A(2 from zero measured using 



standard ESS do not follow the same trend as the those measured through equation 4.2 



This discrepancy suggests that anisotropy is by itself not responsible for the deviations, 
but that they have their origin in measurement noise and error. The mean value of the 
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scaling exponents, excluding the two outliers, was 0.70±0.03. This value is the same as 



the one found by Benzi et al. ( 1995 ) in a wind tunnel at comparable Reynolds number, 
and the value is indicated with a horizontal line in the figure. 



Shen & Warhaft ( 2002 ) found that the exponent of the second order transverse struc- 



ture function was 0.1 smaller in the direction normal to the mean fiow in a wind tunnel, 
both with and without shear. Wind tunnel turbulence is approximately axisymmetric in 
the absence of shear, with the axis of the tunnel corresponding to our z. Although one of 
our data points, namely the one for (Jz/cTr = 1-98, does fall at —0.1, in agreement with 



Shen & Warhaft (2002), their result cannot be reconciled with the body of our data. 

Let us consider the Kolmogorov constant. We calculated the ratio /C2^^ in two 
ways. First, we computed an average of the structure function ratio, {D zzi'f^) / Dr2r2{^)) i 



shown in figure 10 The average was taken over a range of separations that bracketed the 
peak in the compensated structure functions, or 4mm < ri,z < 40mm. In the inertial 



range, according to equations 1 1 . 2 1 and 1 1 . 3[ 

(4.4) 



when n = Q — Q 5 because the dissipation, e, is a scalar quantity. We used 
the standard deviation in the value of the ratio as a measure of the error of the mea- 
surement. The second method was designed to mimic the one usually employed to mea- 
sure C2 (or e) when data are collected in only one direction. That is, we estimated 
each C2 from the maximum of the compensated structure functions. In order to re- 
duce the infiuence of noise, we found the maximum of a polynomial function, second- 
order in log(r), fit to the data between 4 and 40 mm. The ratio between the maxima, 
max(L)^^ (r) /r^/^)/ msD^^Dr^r^ (^) Z^^'^^) , 

was then calculated. As before, the quantity is 
an estimate for the ratio of Kolmogorov constants because the dissipation rate cancels. 
We used the deviation of the data from the fitted parabolas to estimate the error. Al- 



Experimental study of the influence of anisotropy on turbulence 

1.2| ' ' ' ' ' 



25 







1.15 
1.1 
1.05 
1 

0.95 



m 

KH 



0.5 0.6 0.75 1 1.2 1.5 2 

az/(Jr 

Figure 13. The ratio between Kolmogorov constants measured in different directions, 
C^'VC's"^- The O symbols are the ratio between the maximum values of the compensated radial 
and axial structure functions. The □ symbols are {Dzz{ri)/Dr2r2{z)) , the ratio of radial to axial 
structure functions averaged over a range of scales 4 < ri, z < 40 mm. 

though each method described here suffers from limitations, they both serve as practical 
definitions of C^"^^ /C^'l 



13 



The ratios deviated from the 



The two measures of Ci'^^/C^'^ are plotted in fi sure 
isotropic value by up to 4%, except for one outlier whose value was about 8% higher than 
1.0. The observed dependence on the anisotropy, if there was one, was non-monotonic. 



However, Busa et al. (1997) found in a shell model of turbulence that the Kolmogorov 



constant varied monotonically with anisotropy. Although they predicted a dependence 
too weak to detect in our data, it is reasonable to expect one that is qualitatively similar. 
Since we did not see such a monotonic dependence, we concluded that the deviations 
of the Kolmogorov constants from isotropy were probably not due to the anisotropy in 
the fiuctuations, but were rather the result of measurement noise and error. Though 
large, the uncertainty in our data is smaller than the scatter seen in the data collected 



by Sreenivasan (1995), whose amplitude is about 10%. It follows that anisotropy in the 
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velocity fluctuations alone does not account for the variation in previously measured 
values of the Kolmogorov constant. 



5. Conclusions 

We investigated systematically the influence of anisotropic agitation on the inertial 
scales of turbulence in an experiment with Taylor-based Reynolds number Rx = 481. 
Thirty two loudspeaker-driven jets pointed toward the centre of a spherical chamber 
and were driven to produce axisymmetric turbulence, which in a central volume had no 
shear, no mean flow, and was homogeneous. We observed flows for which the ratio of 
axial to radial RMS velocity fluctuations was between 0.6 and 2.3. We found that the 
anisotropy of the velocity fluctuations at the largest scales had the same anisotropy as 
the agitation. According to two inertial range measures, the second order velocity struc- 
ture functions were independent of anisotropy. There was extended self-imilarity, and the 
structure functions had the same scaling exponent in different directions and for different 
anisotropics, 0.70±0.03. The Kolmogorov constant, C2, also showed no dependence on 
the anisotropy. Because of this, we expect anisotropic corrections, such as those uncov- 
ered by the S0(3) decomposition (e.g. Kurien fc Sreenivasan||2QQl[ [Biferale fc Procaccia 



2005), to be smaller than the error in our measurements (about 4%). This holds unless 
the anisotropic corrections cancelled in the two directions we measured, an outcome we 
consider unlikely. Further experiments are necessary to determine whether deviations 
from isotropy present within the uncertainty of our measurement are significant. 
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